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Abstract 

This paper analyzes MIMO systems with multichannel beamforming in Ricean fading. Our results 
apply to a wide class of multichannel systems which transmit on the eigenmodes of the MIMO channel. 
We first present new closed-form expressions for the marginal ordered eigenvalue distributions of complex 
noncentral Wishart matrices. These are used to characterize the statistics of the signal to noise ratio (SNR) 
on each eigenmode. Based on this, we present exact symbol error rate (SER) expressions. We also derive 
closed-form expressions for the diversity order, array gain, and outage probability. We show that the global 
SER performance is dominated by the subchannel corresponding to the minimum channel singular value. 
We also show that, at low outage levels, the outage probability varies inversely with the Ricean if-factor 
for cases where transmission is only on the most dominant subchannel (i.e. a singlechannel beamforming 
system). Numerical results are presented to validate the theoretical analysis. 
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I. Introduction 

Multiple-input multiple-output (MIMO) communication systems have received considerable 
attention in recent years as they offer substantial capacity improvements over conventional single- 
input single-output (SISO) systems [1,2], with no penalty in either power or bandwidth. When 
perfect channel state information (CSI) is available at both the transmitter and receiver, it is 
well known [1] that the capacity-achieving strategy is to transmit on the eigenmodes of the 
MIMO channel using linear transmit-receive processing (hereafter referred to as multichannel 
beamforming (MB)), upon which independent Gaussian codes with water-filling power allocation 
are employed. In practice, the high complexity requirements of Gaussian-like codes are often 
prohibitive, and either suboptimally-coded or uncoded transmission is used. Interestingly, for these 
suboptimal systems, it has been shown that under most performance criteria of practical interest 
(e.g. symbol error rate (SER), mean-square error (MSE), among others), MB still corresponds to 
the optimal choice of linear transmit-receive processing (in some cases up to a rotation matrix) 
[3]. 

In this paper we provide an analytical investigation of the performance of uncoded MB MIMO 
systems in Ricean fading channels. Although these systems are particularly appealing from a 
practical point of view, there are currently very few related analytical performance results available 
in the literature. In [4], the global SER (i.e. SER averaged over all subchannels) was derived for 
Rayleigh and Ricean fading channels under the assumption that all available subchannels were 
used for transmission, with equal power, and with the same (BPSK) modulation. It was shown 
in [5] however (via Monte Carlo simulations), that the performance is significantly improved by 
transmitting on only a subset of the available subchannels, and also by using different powers 
and constellations on each of these. The global SER results of [4] cannot be easily generalized to 
these important scenarios. In [6], global and per- subchannel SER expressions were presented for 
MB MIMO systems in Rayleigh fading in the high signal to noise ratio (SNR) regime, allowing 
for transmission on a set of arbitrarily selected subchannels. In this paper we consider the more 
general class of Ricean fading channels, and seek performance measures for all SNRs. 

The main difficulty in obtaining analytical performance results for MB MIMO systems in Ricean 
fading is that it requires the marginal statistical distributions of the ordered eigenvalues of complex 
noncentral Wishart matrices. Although many results are available for the eigenvalue statistics of 
complex central Wishart matrices (see [7-13], and references therein), there are very few results 
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for the noncentral case. In [14], the joint probability density function (p.d.f.) of these ordered 
eigenvalues was obtained for the special case of full-rank non-centrality matrices. This joint p.d.f. 
was used in [15] to derive the exact p.d.f. and cumulative distribution function (c.d.f.) of the dual- 
antenna MIMO capacity in Ricean channels with rank-1 mean matrices. Unfortunately the joint 
p.d.f. involves a hypergeometric function of matrix arguments and Vandermonde determinants, and 
is not easily marginalized. In [16], the p.d.f. and c.d.f. for the maximum eigenvalue was obtained 
for the special cases of rank-1 and full-rank non-centrality matrices. These results were used 
to analyze the statistics of the output SNR of MIMO maximum-ratio combining (MIMO-MRC) 
systems in Ricean channels. 

In this paper we derive new exact closed-form expressions for the marginal distributions of all 
of the ordered eigenvalues of complex noncentral Wishart matrices. The results apply for non- 
centrality matrices of arbitrary rank. Explicit expressions are given for the marginal c.d.f.s, from 
which the closed-form marginal p.d.f.s can be obtained trivially via differentiation. These marginal 
c.d.f.s are used to analyze the performance of MB MIMO systems in Ricean channels with mean 
matrices of arbitrary rank. New exact expressions are presented for both the subchannel SERs and 
the global SER. These expressions are general in the sense that they allow for transmission on any 
arbitrary number of subchannels, with possibly unequal signal constellations, and with possibly 
unequal powersQ. The exact c.d.f. expressions are also used to obtain new exact closed-form 
expressions for the outage probability of MB MIMO systems. 

In addition to the exact marginal distributions, we also derive new first-order asymptotic ex- 
pansions for the marginal c.d.f.s and p.d.f.s of the ordered eigenvalues of complex noncentral 
Wishart matrices. These expansions are particularly useful for gaining further insights into the 
effect of various system parameters on the performance of MB MIMO systems. In particular, the 
asymptotic p.d.f. expansions are used to derive explicit expressions for the diversity order and 
array gain, which are the factors governing the SER performance at high SNR. These expressions 
reveal that the global SER is dominated by the subchannel corresponding to the minimum channel 
singular value. The asymptotic c.d.f. expansions are used to examine the outage probability of MB 
MIMO at low (practical) outage levels, when only a single subchannel is selected for transmission 
(corresponding to MIMO-MRC transmission). We find that in this outage regime, the outage 
probability becomes independent of the rank of the channel mean, and varies inversely with the 

'Note that while we allow for unequal powers, we assume that they remain fixed; and we calculate the average SER. In other 
words, we are not solving the waterfilling problem, for which the power levels become functions of the eigenvalues and change 
at the fading rate of the channel. 
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Ricean fT-factor. 

II. System Model 
Consider a MIMO system with m transmit and n receive antennas, modeled as 

y = Hs + n, (1) 

where y G C nxl is the discrete-time received vector, s G c mxl j s the transmitted vector, and 
n G C nxl is the noise vector with i.i.d. entries ~ CJ\f(Q, 1). Also, H G C nxm is the Ricean fading 
channel matrix, which is decomposed as follows [17] 

H = eVKU + efl, (2) 

where H is the (arbitrary) rank-L deterministic channel component satisfying tr (HH*) = mn, and 
H is the random (scattered) channel component containing i.i.d. CJ\f(Q, 1) entries. The parameter 
K is the Ricean K-factor, which is the ratio between the energy in H and the average energy in 



H, and e = 1/ \J K + 1 is a power normalization constant. Note that H in © follows a complex 
matrix-variate Gaussian distribution with mean matrix M = y/KeH. and (column) correlation 
matrix e 2 l n . Adopting standard notation from multivariate statistical theory (e.g. see [18, 19]), 
this distribution is denoted 

H ~ CN n ,m (M, e 2 I n ® I m ) . (3) 

Let us now define 

{HHt n ^ m 
(4) 
n > m 

s = min(n, m), and t = max(n, m). With these definitions, W G C sxs follows a complex 
noncentral Wishart distribution, denoted 

W~W a (t,£,n), (5) 

where £ = e 2 I s and 

(6) 

E^MtM n > m 

is the arbitrary-rank non-centrality matrix. 



We consider the class of MB MIMO spatial multiplexing systems considered in [3,5,6]. As in 
[3,5,6], we assume that perfect CSI is known at both the transmitter and receiver. The transmit 
vector can be written as 

s = Bx, (7) 

where B e C mxr is the transmit precoder matrix which maps the r ^ min (to, n) modulated 
data symbols Xj (elements of x, with E {xx^ } = I r , and chosen from possibly different signal 
constellations), onto the to transmit antennas, and is normalized according to 

E {Il s l| 2 } = tr (BB*) ^ P (8) 

where P is the average signal to noise ratio (SNR) per receive antenna. The estimated vector at 
the receiver is given by 

x = Aty, (9) 

where A* e C rx ™ is the receive (spatial) equalizer matrix. 

It was shown in [3], that under many practical design criteria (such as maximizing the mutual 
information, minimizing the arithmetic or geometric mean-square error, among others) the optimal 
transmit and receive filters result in a MB system, and are given by 

B = U H P (10) 

and 

A = (HBB t H t + I n ) 1 HB (1 1) 



respectively, where \J H G C nxr has as columns the eigenvectors corresponding to the r largest 
eigenvalues of H^H, and P = diag (y{^/lh] i=l r j is a power allocation matrix, with ■p i > and 

With this choice of linear transmit and receive filtering, the MIMO channel is decomposed into 
r parallel scalar (eigenmode) subchannels, which are described as follows 



= K k (ey/^kPkXk + ni^ , fc = l,...r (12) 
where n k is a constant (which does not affect the received subchannel SNR), x k and n k are the 
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kth elements of x and n respectively, and 4> k is the kth largest eigenvalue of 



s = £- 1 w ~ w s {t,i s ,n). 



(13) 



The instantaneous received subchannel SNR is given by 



Ik = £ 2 4>k Pk 



k — 1, . . . , r 



(14) 



Clearly the SNR (and therefore the performance) for each subchannel, as well as the overall 
received SNR (and global performance), depend explicitly on the marginal statistical distributions 
of the eigenvalues <pi > ... > <p r of the complex noncentral Wishart matrix in (fl"3l) . In the 
following section we will present new closed-form exact and asymptotic expressions for the 
marginal distributions of these eigenvalues. These results will then be used in Section [IV] to 
analyze the performance of MB MIMO systems in Ricean fading channels. 

III. New Statistical Properties of the Ordered Eigenvalues of Complex 

Noncentral Wishart Matrices 

A. New Exact Ordered Eigenvalue Distribution Results 

In this subsection we derive new exact closed-form marginal eigenvalue c.d.f. expressions. 
Note that exact marginal p.d.f. expressions can also be easily obtained from these c.d.f. results 
via differentiation. These results, however, are omitted due to space constraints. For convenience, 
we consider the smallest, largest, and kth. largest eigenvalues separately. The proofs of all results 
in this section are given in the appendices. 

First consider the smallest eigenvalue <j) s . It should be noted that, in additional to the performance 
analysis of MB MIMO systems considered in this paper, the statistical properties of the smallest 
eigenvalue of Wishart matrices are important in the analysis of various other MIMO systems and 
applications (see e.g. [20-23]). 

Theorem 1: The c.d.f. of the smallest eigenvalue S of the complex noncentral Wishart matrix 
S in (fT3l is given by 



f^ s (x) = i- |*(*)|/ 1* (o)|, 



(15) 



where ^ (x) is an s x s matrix function of x 6 (0, oo) whose entries are given by 




(16) 



and where Ai > . . . > Ax are the L s) non-zero eigenvalues of Vt. Also, Q Ptq (a,b) is the 
Nuttall Q-function, defined in [24] by 

f°° ( x 2 + a 2 \ 
Q P ,g (a, b) = J x p exp ( J I q (ax) dx, (17) 

I q (■) is the gth order modified Bessel function of the first kind, and T (•, •) is the upper incomplete 
gamma function, defined as [25] 

/co P^ 1 k 

t^e-Mt = (p - 1)! e- J2 Tf. p = 1, 2, . . . . (18) 

Proof: See Appendix [I] □ 

Note that the result in (fT5l) can be easily programmed and efficiently evaluated. Moreover, the 

sum of the Nuttall (^-function indices in (fT6l) is odd. As such, this function has a closed-form 

representation given in [26, Eq. 8]. 

Now consider the largest eigenvalue fa of S. In [16], the c.d.f. of fa was derived for the 

particular cases where the mean matrix was either rank-1 or full-rank, and used to analyze the 

performance of MIMO-MRC systems in Ricean channels. The following theorem presents a new, 

simpler expression for this c.d.f., and also generalizes the results of [16] as it applies for mean 

matrices with arbitrary rank. 

Theorem 2: The c.d.f. of the largest eigenvalue fa of the complex noncentral Wishart matrix 

S in (fT3l) is given by 

F 01 (x) = |H(x)|/|*(O)|, (19) 
where S (x) is an s x s matrix function of x G (0, oo) whose entries are given by 



{s(s)K 



■J 



2<*—t)/2 [Q s+t _ 2l+u _ s (y2A~ 0) - Q s+t -2i+i,t- s ( V2x)} J = 1, • • • , L 

j(t + s-i- j + l,x) j = L + l,..., 

(20) 

where 7(-, ■) is the lower incomplete gamma function, given by 

i-x ( P _1 k\ 

7 (p,x) = j t^ 1 e- t dt=(p-l)\ h-e-^^j, p=l,2,.... (21) 
Proof: See Appendix HB □ 
Finally we consider the A;th largest eigenvalue <p k of S. 

Theorem 3: The c.d.f. of the A;th largest eigenvalue cj) k of the complex noncentral Wishart matrix 
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S in (fT3lh where k > 1, is given by 

F^ (x) = (x) + Pr ((f> 8 < . . . < x < < . . . < (j>x) = F 0fc _ 1 (a;) + p, (22) 



where 



p = c 3 J2\ @ ( x )\> ( 23 ) 



(2L-s-t)/2 



c 3 = r , (24) 

f {*(x)V , z = l,...,Jfe- 1 

{ (*)U,-= (25) 

indicates the summation over the combination (ai < a 2 < ■ ■ ■ < a fc _i) and (a fe < a k+i < . . . < a s ), 

i 

(ai, • • • , «s) being a permutation of (1, . . . , s). 

Proof: See Appendix Hill □ 
Remark: Let u>i ^ cu 2 ^ • ■ ■ ^ c<j s ^ be the ordered singular values of H. Recalling that 
S = E _1 W, we have the following relationship 

F u . (x) = (e- 2 x 2 ) . (26) 

B. New Asymptotic Ordered Eigenvalue Distribution Results 

In this subsection we present new asymptotic first-order expansions of the marginal eigenvalue 
p.d.f.s and c.d.f.s. of complex noncentral Wishart matrices. These will be particularly useful for 
deriving the diversity order and array gain of MB MIMO systems in the following section, as 
well as for analyzing the asymptotic outage probability. Note that a first order expansion of the 
p.d.f. of the A;th largest eigenvalue of complex central Wishart matrices was presented in [6]. 

Theorem 4: The first order expansions of the marginal p.d.f. and c.d.f. of the /cth largest 
eigenvalue 4> k of the complex noncentral Wishart matrix S in (fT3l) . where 1 < k < s, are given 
respectively by 

U k M = a k <j> d k k + o (<Pt k ) ( 27 ) 

and 



r^<>u=^^ +1 + o(<Pt +1 ). (28) 
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where 



dk — s k t k — 1 



(29) 



and 



Sfc = s — k + 1, tk = t — k + 1. 



(30) 



Also, a*,, is given for = 1 by 



T s (s) 

r s {t + s ) 



-tr(fi) 



(31) 



and for > 1 by 



s k t k Tk-i(s)T Sk (s k ) 

a k = j n c 3 

T Sk (t k + S k ) 



(it, 



A 



(•-0/2 



IXI 



(32) 



where L.(-) is the normalized complex multivariate gamma function, defined in (|48T) . and X is an 
s x s matrix with (2, j)th entry 



{X} 





i = 


1,. 




1. 


j = l,...,L 






i = 


k,. 


.,s, 


j 


= 1,...,L 




r-i) 


i = 


1,- 


.,k- 


1, 


j = L + l,... 


,s, j >i 





i = 


1,- 


.,k- 


1, 


j = L + l,... 


,s, j <i 


(_l)i-j( a _j)! 


i = 


fc,. 


.,s, 




= L + l,...,s, 


J < i 





i = 


fc,. 


.,s, 


J 


= L + l,...,s, 


J > i 



where 



4"> M 



i=0 



+ (— x) 1 
k — i 



is the generalized /cth-order Laguerre polynomial. 
Proof: See Appendix [IV] 



(33) 



(34) 



□ 



IV. Performance Analysis of MIMO Multichannel Beamforming in Ricean 

Fading Channels 

A. Symbol Error Rate Analysis 

We now analyze the SER performance of the MB MIMO systems introduced in Section HB 
For many general modulation formats (see below), the average SER of the A;th subchannel can be 
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expressed as [27] 

SER fc = E lk {a fe Q (VWklk) } , k=l,...,r (35) 

where Q(-) is the Gaussian Q-function, and a k and (3 k are modulation- specific constants. Some 
example modulation formats for which (|331) apply include BPSK (a k = 1, = 1); BFSK with 
orthogonal signalling (a k = 1, (3 k = 0.5) or minimum correlation (a k = I, (3 k = 0.715); and 
M-ary PAM (a k = 2(M - 1)/M, (3 k = 3/(M 2 - 1)). Our results also provide the approximate 
SER for those other formats for which (|35l) is an approximation, e.g. M-ary PSK (a k = 2, (3 k = 
sin 2 (7r/M)) [27, Eq. 5.2-61]. Using results from [10] and [28], (1351) can be expressed in the 
following equivalent form 

SER fc = / — -5— dw, fc = 1, . . . , r (36) 

2V^ Jo WPkJ 

where we have used the fact that F Jk (u) = (e~ 2 u/p k ). Hence, by applying Theorems \1§3\ in 
(l36l) we obtain an exact expression for the average SER of each subchannel. Although it does not 
appear that the integrals in d36l) can be evaluated in closed form, numerical integration can be 
performed to evaluate SER fc much more efficiently than is possible via Monte Carlo simulation. 

The global SER can be derived from the subchannel SERs as follows; since independent symbols 
are sent on each subchannel during each channel use: 

1 r 

SER = - 2j SER fe . (37) 

r k=l 

To gain further insights, we now consider the SER at high SNR. We will restrict this asymptotic 
analysis to systems with uniform power allocation (i.e. p k = P/r), since (as mentioned in [6]) 
most of the practical power allocation solutions in [3] tend to uniform as the power is increased. 
In the high SNR regime, the key factors governing system performance are the diversity order 
and array gain. We will now present closed-form expressions for these factors. 

Using a general SISO SER result from [29], we find that in our case SER^ can be approximated 
in the high SNR regime by considering a first order expansion of the p.d.f. of <\> k as <\> k — > + . 
Hence, using the result from [29], along with Theorem,^ we obtain a high SNR subchannel SER 
expression given by 

SER^ = (G a (k) ■ py G ^ + o (P-^W) , k = 1, . . . , r (38) 
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where the diversity order is 



G d (k) = 4 + 1 



(39) 



and the array gain is 



2f3 k e 2 f a k 2 d «a k r(d k + 3/2) 
r V yft(d k + l) 



) 



i/(d fc +i) 



G a (k) 



(40) 



Comparing (l39l with the i.i.d. Rayleigh results presented previously in [6], we see that the 
subchannel diversity orders are the same in both Rayleigh and Ricean channels. Moreover, since 
Grf(l) ^ Gd(2) ^ ••• ^ Gd(r), the rth subchannel has the poorest performance in terms of 
average SER. Using (|37|). the global average SER of MB MIMO systems at high SNR can be 
obtained as 



which is clearly dominated by the rth subchannel SER (i.e. the subchannel corresponding to the 
smallest singular value). 

B. Outage Probability Analysis 

We now consider the outage probability of MB MIMO systems in Ricean fading channels. 
The outage probability is an important quality of service measure, defined as the probability that 
the received SNR drops below an acceptable SNR threshold 7th- For convenience, we assume 
an equal power allocation strategy is employed. In this case, the subchannel SNRs are ordered 
according to 71 > . . . > 7 r in (fT4l) . and the outage probability of the overall MB MIMO system is 
dominated by the subchannel corresponding to k = r. As such, the outage probability is obtained 
exactly from Theorems 1-3 as follows 



We note that for the special case r = 1 (i.e. only the best subchannel corresponding to 0i is 
used), the MB MIMO system we consider is equivalent to the MIMO-MRC systems considered 
in [11, 16]. For these systems, outage probability expressions were derived previously in [16] for 
the special case of channels with rank-1 and full-rank mean matrices. Our result (|42l) is clearly 
more general as it applies for all r ^ 1 and for mean matrices of arbitrary rank. Moreover, for 
the special case r = 1, our result is simpler than the previous results given in [16]. 



SER°° = - (G a (r) ■ P) 



G d (r) 



+ o(P- g ^) 



(41) 




(42) 
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In practice, we are usually interested in small outage probabilities (i.e. 0.01, 0.001, ...), which 
correspond to small values of 7^. To gain further intuition at these small outage probabilities, let 
us consider the case r = 1 (i.e. the MIMO-MRC case), and use ([28]) in Theorem |?] to write the 
outage probability in (1421) as follows 

~ T s (s)e- tI ^ f lth (K + iy 



(^^)%o(( 7th n. (43) 



° UI T s (t + s 

Since tr (HH*) = st (see Section HU) we can further simplify to obtain 



Pout = Tjt + s) Vp) + ° ((7th) } • (44) 

This explicitly shows that for MIMO-MRC transmission the outage probability (at low outage 
levels) does not depend on the rank of the channel mean. Moreover, since 

4^ (K + l) st e~ Kst = -Kste- Kst (K + l) 5 *" 1 < 0, K > (45) 
dK 

we see that the outage probability varies inversely with the Ricean if-factor. 

V. Numerical Results 

For our numerical results we consider a 3 x 5 Ricean MIMO channel and, unless otherwise 
specified, a rank-3 deterministic component H with singular values {2.9751, 2.2840, 0.9657}, 
which were randomly generated to verify the analysis. 

It is important to note, however, that all of the analytic results in this paper apply for arbitrary 
antenna configurations, and for arbitrary deterministic channel components. 

Fig.[TJgives the marginal c.d.f.s of the ordered channel singular values. The analytical curves are 
generated using Theorems 1-3 and the singular value relationship (|26l) . and the simulated curves 
are generated based on 100,000 channel realizations. The figure shows perfect agreement between 
the analytical results and the simulations. 

Fig. |2] shows the c.d.f. of the smallest singular value for different Ricean /^-factors, and for 
rank-1 and rank-3 mean matrices. The deterministic component for the rank-1 case was generated 
based on the channel model in [15]. As expected, for both mean matrices, as K becomes small 
(i.e. K = — lOdB), the c.d.f.s converge to that of the smallest singular value of a Rayleigh channel. 

Fig. |3] shows the exact subchannel SER curves based on (|36l ), exact global SER curve based 
on d37]), and Monte-Carlo SER simulation results, for a MB MIMO system with K = OdB. All 
subchannels are used with BPSK modulation (a^ = 1, (3k = 1) and uniform power allocation. 
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High SNR curves based on (|38l ) and (I4TT) are also presented. In all cases, there is exact agreement 
between the analytical SER results and the Monte-Carlo simulations, and the diversity order and 
array gains predicted by the high SNR analytical results are accurate. We also see that the SERs 
of the 1st and 2nd subchannels are significantly better than the 3rd subchannel SER. This suggests 
that further performance improvements may be possible by using only a subset of the subchannels 
for transmission (with higher order constellations). The following figure investigates this further. 

Fig. |4] shows the analytical and Monte-Carlo simulated global SER curves for MB MIMO 
systems with different numbers of r active subchannels. The analytical curves are generated based 
on (1371) . Uniform power allocation is assumed and, for a fair comparison, the overall rate is set 
to 3 bits/s/Hz in each case. For r = 1, 8PSK («x = 2, Pi = 0.146) is employed; for r = 2 we use 
QPSK {pii — 2, Pi — 0.5) for the first subchannel and BPSK for the second subchannel; and for 
r = 3 we use BPSK for all subchannels. For the r = 3 case see an exact agreement between the 
analytical and Monte-Carlo simulated curves. As discussed in Section IIV-AI (1371 ) only provides 
an approximation for QPSK and 8PSK, however we see from the r = 2 and r = 3 curves that the 
approximation is accurate. In particular, for low to moderate SERs (i.e. SER < 10~ 3 ) the analytical 
curves match almost exactly with the simulated curves. We also see that the SER for the r = 1 
and r = 2 cases is significantly better than for the case where all subchannels are used, which 
is in agreement with previous Rayleigh fading observations in [5] (via Monte-Carlo simulations). 
Moreover, we see that the r = 1 system has a higher diversity order than the r = 2 system (since 
we've shown the diversity order to be dominated by lowest singular value subchannel), but is 
shifted to the right due to the higher order constellations. This motivates the design of practical 
adaptive subchannel selection algorithms, which is an interesting topic for future work, but beyond 
the scope of this paper. 

Fig. [5] shows analytical and Monte-Carlo simulated outage probability curves for a MB MIMO 
systems with r = 1 (i.e. MIMO-MRC transmission), comparing different Ricean ^-factors. The 
analytical curves are generated based on (1421) . We see an exact agreement between the analytical 
and simulated curves in all cases. We also see that increasing the Ricean if-factor results in a 
reduction in outage probability (and an improvement in system performance) at low outage levels. 
This agrees with the analytic conclusions given in Section IIV-BI It is also interesting to observe 
that the opposite occurs in the high outage regime (i.e. increasing the fC-factor increases the 
outage probability). 
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VI. Conclusion 

We have examined the performance of MEMO systems employing multichannel beamforming 
in arbitrary-rank Ricean channels. Our results are based on new closed-form exact and asymptotic 
expressions which we have derived for the marginal ordered eigenvalue distributions of complex 
noncentral Wishart matrices. We have presented exact and high-SNR SER expressions, and derived 
the diversity order and array gain. Our results have shown that the global SER performance is 
dominated by the subchannel SER corresponding to the minimum channel singular value. We 
have also derived new closed-form exact expressions for the outage probability and, for the case 
of MIMO-MRC transmission, have shown that for outage levels of practical interest, the outage 
probability varies inversely with the Ricean K-factor. 

Appendix I 
Proof of Theorem Q] 

We require the following Lemma, which gives the joint eigenvalue p.d.f. of S for the case 
where the non-centrality matrix O has arbitrary-rank. 

Lemma 1: The joint p.d.f. of the ordered eigenvalues (j>\ > </> 2 > . . . > 4> s > of the complex 
noncentral Wishart matrix S in (fT3l is given by 

s s 

f . . . , <j>.) = d|T| fj {<f>i - <Pj) J] <Pl- s e-^, (46) 

i<j k=l 



where 



and where 



c = ""^ {{t - (47) 

' 1 r s _ L (s - l) nil (\r L ) nL (a, - a,) 



r s (t) = H(t-i)\. (48) 



1=1 



Also T is an s x s matrix with (i,j)th entry 



m..= I v . J ' (49) 



3 F X (t - s + 1; Xj<f>i) i = l,...,s, j = l,...,L 

(t-jy- 

where o-^i(') is the scalar Bessel-type hypergeometric function. 

Proof: We start by combining a result from [14], which gave the joint eigenvalue p.d.f. for 
the special case of full-rank Q, with a hypergeometric function determinant result from [30], to 
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express the joint eigenvalue p.d.f. in the full-rank ft case as follows 



FR {<Pl,...,<p s 



p-tr(n) p -tr(*) » 

)= — —. n 



■t-8 10 



F l (t-s + l;X 



(50) 



((t _,,,). 11- "'ii- n*<,(A,-A J ; 

where A = diag (Ai, . . . , X s ), Ai > ... > X s > are the eigenvalues of f2, and $ = diag (<f>x, . . . , S ). 
We generalize this result to arbitrary-rank matrices fi by taking limits 

e -tr(n) e -tr(*) 



/(0i,...,0 a ) = lim / F r(0i, ••■ 

-\l-i-i,...,A s — >0 



where 



lim 



((t- s )\y 



\fi(*j)\ 



n 



U^C (51) 



fc=i 



(52) 



with fi (Xj 



A i+ i,...,A,-0 n'<i (Ai - Ay) 

o-Fi (i — s + 1; Xj4>i). To evaluate these limits we apply [31, Lemma 2] to obtain 

MAO ... fi (Ai) /{-^(o) ... A (0) (0) 



lim 



l/i(Ai)| 



/• (Ai) ... f s (X L ) f { s M) (0) ... /i 0) (0) 



a, +1 ,..,a^o n:<, (Ai - \j) n ^ {Xi _ A .) ^ A r L )r s _ L (* - l) 



(53) 



where the required derivatives are easily evaluated as 



if (o) 



(t-s + k)\ 

The result now follows by substituting (I52l)-(l54l) into (I5TT) and simplifying. 

We now proceed to evaluate the c.d.f. of the minimum eigenvalue 4> s as follows 



(54) 



□ 



F 0s (x) = 1-Pr(0x > ...0 S >ar) = 1- / /(&,...,&) d& -d& (55) 
where £>! = {a; < S < . . . < </>i}. To evaluate the integrals in (1551) we require the following result 

i*i n - 0i) = i*i kr] = E E s § n o*) fi c fe {*u^ (56) 



k=l 



which is easily obtained using Lemma [7] and the definition of the determinant. In (1561) , /x = 
(/xi, /x 2 , . . . , /x s ) and a = (<7i, cr 2 , • • • , cr s ) are permutations of (1, . . . , s), the sums are over all 
permutations, and sgn(-) denotes the permutation sign. Substituting (|46|) and (l56l) into (1551) . and 
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using [12, Lemma 1], it can be shown that the c.d.f. of <p s can be written as 



{x) = 1 - Pr (<f> s > x) = 1 - ci 



fi(Xj,y) dy 



(57) 



where 



Using the relation 



y^e-foFi^-s + ljAj-j/) j = l,...,L 

yt+ s-i-j e -y (j- - s )\/(t j — L + 1, . . . ,S 



(58) 



jFi (* - s + 1; x) = (t - s)!x-^ s )/ 2 /^ s 



as well as (TT7T ), the remaining integral in (1571) can be evaluated as follows 



(59) 



j = L + l,...,s 



(60) 



(t-s)!r(t + s -2-j + l,x)/(t-j)! 
Since Pr (</> s ^ 0) = 1, Ci can also be expressed as 

ci = 1/|Xm(0)|. (61) 

Substituting (|60l) and (|6TT) into (1571 ) and simplifying by removing common factors from the 
numerator and denominator determinants, we obtain the desired c.d.f. of the smallest eigenvalue. 

Appendix II 
Proof of Theorem [2] 

We can evaluate the c.d.f. of the maximum eigenvalue <\>\ as follows 



Ffr (x) = Pr {<f> 3 < . . . < X < x) = / / 

Jv 2 



> i . . . . , (p s ) d(j)i 



(62) 



where D 2 = {<p s < . . . < 4>\ < x}. We evaluate these integrals by following a similar procedure 
to that used for evaluating (l55l) in the proof of Theorem [Q to obtain 



Ffa (x) = c 2 



fi (Xj,v) dy 



(63) 
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where (Aj, £/) is defined as in (|58l) . We obtain the desired result by noting that C\ = c 2 since 
Pr (4>i < oo) = 1, and using the property 

X r+oo f+oo 

f i (\ j ,y)dy= / f i (\ j ,y)dy- / f i (X j ,y)dy. (64) 

J0 Jx 



along with d60|) and fCD in (1631) . 

Appendix III 
Proof of Theorem [3] 

Let D 3 = {4> s < . . . < x < (f) k -i < . . . < 4>\}, D 4 = {x < 4> k _i < . . . < (j)i < +00} and D 5 
{0 < 4> s < . . . <pk < x }- Using (|46l ) and (|56l ) we can write the probability p as 



p = / e e s § n ^ n ^ fce ~^ m^^- < 65) 

Da cr /i fc=l 

Note that the summation over a can be decomposed as [12] 

£=E£, £ r w 

cr 1 

where denotes summation over the permutations (r ai , . . . ,r ak l ^ of (1, . . . , k — 1) and 

X] rQ2 denotes summation over the permutations (r afc , . . . , r a J of (A;, . . . , s). Therefore, we have 

1 M Ds fc=l 

1 /1 

where 

, re 1 re— 1 /»-)-oo 



Ii fa, 

D 4 



/ p+00 
n m ra< ^ = n / cr^- ro^ ^ - «*) 
,, i=i i=i j:i ' 



D% — k % — k 



The last equality follows from [12, Lemma 1]. The desired result follows from (1641) . (1601) and the 
definition of the determinant. 
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Appendix IV 
Proof of Theorem U 

Here we will derive the p.d.f. expansion (1271) . The corresponding c.d.f. expansion (|28l) then 
follows trivially. 

We start by noting that since (<p k ) = dF^ (<f) k ) /d(j) k , the Taylor expansion of f^ k (<p k ) 
around the origin can be written as 



p k M = f£ (o) + (o) <p k + ■ ■ ■ + 



(70) 



In order to simplify the derivations, we will initially work with F^ix) under the assumption 
of full-rank Q (i.e. L = s). We will then generalize our result to the arbitrary-rank Ct case by 
evaluating limits where necessary. 



A. Derivation for <pi 

We first derive the first order expansion of ((pi). Using ([T9l) and a well-known result for the 
A;th derivative of a determinant, we obtain 



F[ q+l) (x) 



1 



E 



(g + l)! 



d 9 * 



dx qi 



{3(*)} 



'■j 



(71) 



91 H Vq s =q+l 

We require the gjth order derivative of the Nuttall Q-function Q s +t-2i+\,t-s (y/2Xj, V2x) in (|20]>. 
Omitting details, with the help of Leibnitz' rule, we evaluate these derivatives as follows 



r=0 



(a:) (x) 



for qi ^ 1, where 



(x) = exp (-Xj - x) 



92 {x) 



A 'J l V2 J 2-, p! (t_ s + p )! 



(72) 

(73) 
(74) 



Hence, 



(«) 

s+l— 2i+l,t— a 



x=0 







r=t—i 



(r-t+i)\(r-s+i)\ 



qi — 1 < t — i 



qi — 1 ^ i — i 



(75) 
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To obtain the first order expansion of (0i), we require the minimum exponent q in (1701 ) 
(and corresponding values of qi,...,q s in (1711)), such that (0) 7^ 0. Using dTTI ), (1751 ) and 

the properties of the determinant, we find that 



where i — 1, 



qi-1 = t- i+pi, 
s and (pi, . . . ,p s ) is a permutation of (0, . . . , s — 1). Thus 

s 

d\ = q = qi — 1 = st — 1 , 



(76) 



and 



where 



{Ax}-, , = i/(t + s-i-j + l). 



(77) 



(78) 



(79) 



We now simplify a\. Using [26, Eq. 13] we can write 



L 



(t-s) 



-A,) 



(80) 



i*(o)i=2-(*- 1 )/ 2 ri A ^ )/2 

\i=0 / \j=l 

We also manipulate |Ai| by subtracting the first row from all other rows, removing factors, then 
subtracting the first column from all other column, and again removing factors. This yields 

N 2 



|Ai| 



((s - 1)\ (t - iy.y [ lfl] 



(81) 



(t + s-l)\(t + s-2)V 

where A J 1 ' 1 ' is the principle submatrix of Ai, with first row and column removed. Continuing 
this same process s — 1 more times, we obtain 



|Ai| 



n^ + s-i)' r i+s (t + s ) 

Substituting ([80]) and ([82]) into d78]) and simplifying yields 

ste-*^ Yl^. (A, - A,) _ ste^ | A /-i| 



(82) 



ox 



r s (t + s) 



,(t + s) \Lfjf\-\j)\ 



(83) 



Finally, we remove the remaining determinant ratio. To do this, we manipulate the numerator 
determinant in such a way as to construct a scaled version of the denominator determinant. Start 
by considering the row % — s. Using (l34l) . the elements in this row of the denominator determinant 
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can be written as 

s - 2 ' I 1 \A> A, 1 *-" 



L 



(t-s), 
s-1 \ 

k=0 



**> = E Ui-Jl+CT '' = 1 -- (84) 



We construct this row as the last row of the numerator determinant by first dividing the sth row in 
the numerator determinant by (s — 1)!, and multiplying the determinant by (s — 1)! to compensate. 
This gives elements corresponding to the (s — l) th order polynomial terms on the right-hand side 
of (1841) . All of the remaining terms on the right-hand side of ([841) can be generated from weighted 
sums of the rows i = 1, . . . , s — 1 in the numerator determinant (these row operations do not 
change the value of the determinant). If we then apply the same process to rows i = s — 1, . . . , 1, 
in order, at each stage constructing another row of the denominator determinant, and pulling out 
a multiplicative factor of (i — 1)!, we obtain 



i \ * — 1 1 s 



Y[(s-i)\ = T 8 (s). (85) 



Substituting (|85l) into (|83l) gives the desired result for a\. Note that in this case, a\ is only a 
function of tr(fi), and is independent of the rank of fi. 

B. Derivation for 0*. (k = 2, . . . , s) 

Now consider the case k > 1. Starting with (1221) for the case of full-rank ft, We observe that 
the minimum exponent q is obtained when (ax, ... , a s ) = (1, . . . , s). Applying the same steps as 
for the (pi case, we obtain 




(86) 

...,8 

where (e^, . . . , e s ) is a permutation of (0, . . . s — k). Hence 

s 

4 = q - 1 = 7^ff» - 1 = s k t k - 1, 2 ^ k < s. (87) 

i=k 



Combining ([77]) and ([87]) yields (1291 . Also, we obtain 

r s fe (^) r Sfe (^. 



Sfc t fc c 3 2^^- 1 )/ 2 |A 2 ||A 3 | 
= — ? 71 \ (88) 



where 



Q s+t-2i+i,t-s (a/2Aj, 0) i — 1, . . . , k — 1 
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and 

{A 3 } id = l/(t k + 3 k -i-j + l) i,j = l,...,8-k + l. (90) 



Now, using [26, Eq. 13], we factorize |A 2 | as follows 



where 



We evaluate |A 3 | using ^T9]) and ((82) as 



\U. X i (t ~ s)/ ^J r fe _ 1 ( s )2(- fe2+fe (* +s+1 )-(* +s »/ 2 |A 4 | (91) 



E/-n( « = l,...,Jfe-l 

{a 4 }. . = <; z ^-° ls ;^ j/ (92) 

A! l e~ x i i = k,...,s 



i A i _ (r Sfc (sfc)r Sfc (t fe )) r f)fe _ Sfc (t fc - s fe ) 

I ^3 1 — 77 — : x • yyj) 

1 t k +s k {tk + s k ) 



Now substituting (1931 ) and (1911) into (1881) and simplifying yields 

= gfc*fcr fc _i(s)r Bfc (s fc ) |A 4 | 
r Sfe (4 + Sfc ) rii<i(^- A i) 

Finally, we generalize this result to arbitrary-rank fl. To this end, we require the following limit 

|A 4 | 

l im TTs — 71 rri (95) 

A i+1 ^0,...,A s ^0 Y[ i<:j (Aj - Aj) 



which we easily evaluate using (1531) . and substitute into (1941 ) to obtain the final expression 
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x 



Fig. 1. Ordered c.d.f.s of the singular values of a 3 x 5 Ricean MIMO channel, with K = lOdB. 




Rayleigh 
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x 

Fig. 2. Smallest singular value c.d.f. for a 3 x 5 Ricean MIMO channel with rank-1 and rank-3 
channel means, and for different fT-factors. Rayleigh c.d.f. presented for comparison. 



24 




SNR(dB) 

Fig. 3. Exact analytical SER, high SNR analytical SER, and Monte-Carlo simulated SER for 
3x5 MB MIMO in uncorrelated Ricean fading, with rank-3 mean matrix and K = OdB. 
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Fig. 5. Outage probability of 3 x 5 MIMO-MRC in Ricean channels with rank-3 mean matrix, 



various ^-factors, and for P = OdB. 



